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XMM-Newton discovery of transient 285.4 s X-ray pulsar 
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1 INTRODUCTION 



Since their discovery (|Giacconi et al.lll97ll ). accreting X-ray 
pulsars have been prime targets for both observational and 
theoretical study dwhite. Swank fe HoltJ[l983l ; lNagasdfl989l ; 
iBildsten et al.ll 19971 ). X-ray pulsars proved to be unique lab- 
oratories that allow to test a wide variety of physical pro- 
cesses in the presence of magnetic fields and strong gravity, 
and study stellar evolution in binary systems. 

Until recently, the study of pulsating X-ray sources 
was limited to our Galaxy and the nearby Magellanic 
Clouds. The advanced capabilities of Chandra and XMM- 
Newton observatories opened the possibility to study spec- 
tral and timing properties of individual brig ht X-ray sources 
associated with more distant galaxies (|Fabbianol 120061 ; 
iFabbiano fe Wh ite 2006), allowing to extend searches for 
pulsating X-ray binaries beyond the Milky Way and its 
immediate neighbors. The nearby normal galaxies with 
recent star formation are especially suitable candidates 
for such studies, since they provide the environments ex- 
pected to harbor a significant population of accreting bi- 
nary systems with highly magnetized neutron star and white 
dwarf primaries, some of them bright enough to be de- 
tected by current X-ray instruments. Recent analysis of 
the monitoring observations with XMM-Newton and Chan- 
dra has already led to the discovery of coherent pulsa- 
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ABSTRACT 

We report on the discovery and analysis of the transient X-ray pulsar XMMU 
J013359. 5+303634 detected in the 2010 July- August deep XMM-Newton observations 
of the nearby spiral galaxy M33. The X-ray source exhibits pulsations with a period 
P ~ 285.4 s and pulsed fraction ~47 per cent in the 0.3-10 keV energy band. The 
X-ray spectrum of XMMU J013359. 5+303634 is hard and is well fitted with an ab- 
sorbed simple power law of photon index T ~ 1.2 in the 0.3-10 keV energy band. The 
X-ray properties of the source and the absence of an optical counterpart brighter than 
20 mag suggest that XMMU J013359. 5+303634 is not a Galactic foreground object 
and probably belongs to M33. The maximum observed luminosity of the source in the 
0.3-10 keV energy range is ~ 1.4 x 10 37 ergs s _1 at 817 kpc, at least 20 times higher 
than quiescent luminosity. Based on the transient behavior, overall X-ray properties 
and possible optical counterpart, we classify XMMU J013359. 5+303634 as a Be/X-ray 
binary candidate. 

Key words: galaxies: individual: M33 - X-rays: binaries. 



tions in the t hree X-ray sources located in the Local Group 
galaxy M31 (lOsborne et al.ll200ll; ITrudolvubov et alj 12009 ; 
iTrudolvubov fe Priedhorskvl 120081 ) and two X -ray pulsars 
in the nearby spiral galaxies NGC 2403 a n d NGC 1313 
jTrudolvubov. Priedhorskv fe Cordoval 120071 ; iTrudolvubovl 
2008). These results clearly demonstrate that bright (Lx > 
10 37 ergs s _1 ) pulsating X-ray sources can be detected in 
the galaxies inside and beyond the Local Group up to the 
distances of a few Mpc. 



At the distance of 817 kpc dFreedman et al.ll20oil ). the 
nearest late-type face-on spiral galaxy M33 provides an ex- 
cellent opportunity to study X-ray sources associated with 
young stellar population. M33 was a target of multiple X - 
ray obse rvations with the Einst ein (iTrinchieri et al.lll988l). 
ROSAT dHaberl fe PietscbfeOOll) . Chandra |Plucinskv et all 
l2008t ; iTullmann et alj|201ll), and XMM-Newto n observato- 
ries l|Pietsch et all 12004 iMisanovic et~ai1l2006l ). These ob- 
servations uncovered a population of hundreds of X-ray 
sources with many of them identified as X-ray binary and 
supernova remnant candidates associated with regions of re- 
cent star formation. 

In this Letter, we report on the discovery of the coher- 
ent 285.4 s pulsations in the flux of transient X-ray source 
XMMU J013359. 5+303634 in M33, using archival data of 
XMM-Newton observations. We study X-ray properties of 
the source, search for optical counterpart and discuss its na- 
ture. 
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Table 1. XMM-Newton observations of M33 used in the analysis of XMMU J013359. 5+303634. 



Date, UT 


Obs. ID 


Instrument 


Mode/ 


RA (J2000) a 


Dec (J2000) a 


Exp. 6 








Filter 


(h:m:s) 


(d:m:s) 


(ks) 


2012 Jul. 9-10 


0650510101 


EPIC 


Full/Medium 


01:34:13.40 


+30:47:48.8 


99.7 


2012 Jul. 11-12 


0650510201 


EPIC 


Full/Medium 


01:33:45.22 


+30:35:44.4 


99.6 


2012 Aug. 12-13 


0650510601 


EPIC 


Full/Medium 


01:34:21.29 


+30:27:07.3 


100.0 



a — pointing coordinates 
b - instrument exposure used in the analysis 



2 OBSERVATIONS AND DATA REDUCTION 

To study timing and spectral characteristics of XMMU 
J013359.5+303634, we used the data of 2010 July 9-10 
and July 11-12 XMM-Newton observations of M33 fields 
with three European Photon Imaging Camera (EPIC) in- 
struments (MOS1, MOS2 and pn) ^Turner et all l200ll : 
IStrueder etafll200ll ) (Table [J). We also used several Chan- 
dra (Obs. IDs: 786, 1 730, 6376, 6377 6385 6 387, 6388, 6389 
7170, 7171, 7196) jPlucinskv et all l200Sl : iTullmann etaU 
120111 ) observations with the field of view covering the source 
position to obtain upper limits on the quiescent source flux. 

We reduced XMM data using XMM-Newton Science 
Analysis System (SAS v ll.O.ojB Before generating X-ray 
images, source spectra and light curves, we performed stan- 
dard screening of the original event files to exclude time in- 
tervals with high background levels. The standard SAS tool 
barycen was used to perform barycentric correction on the 
original EPIC event files used for timing analysis. 

We generated EPIC-pn and MOS images of the M33 
fields in the 0.3-10.0 keV energy band, and used the SAS 
standard maximum likelihood (ML) source detection script 
edetect_chain to detect and localize point sources. We used 
bright X-ray sources with known counterparts from USNO- 
B an d LOGS catalogs (|Monet et alj 120031 ; iMassev et ail 
2006) and Chandra source lists to correct EPIC image as- 
trometry. After correction, we estimate residual systematic 
error in the source positions to be of the order ~ 1 arcsec. 

To extract EPIC-pn source light curves and spectra dur- 
ing the XMM-Newton observations, we used the elliptical 
regions including at least ~ 70 per cent of the source energy 
flux. The adjacent source-free regions were used to extract 
background spectra and light curves. The source and back- 
ground spectra were then renormalized by ratio of the re- 
spective detector areas. For spectral analysis, we used data 
in the 0.3— 10 keV energy band. In this analysis we use valid 
pn events with pattern 0-4 (single and double) and pattern 
0-12 (single-quadruple) events for MOS cameras. To syn- 
chronize source and background light curves from individual 
EPIC detectors, we used the identical time filtering criteria 
based on Mission R elative Time (MRT), following the pro- 
cedure described in lBarnard et all (|2007t ). The background 
light curves were not subtracted from the source light curves, 
but were used later to estimate the background contribution 
in the calculation of the source pulsed fractions. 

The EPIC spectra were grouped to contain a minimum 
of 20 counts per spectral bin in order to allow \ 2 statis- 
tics, and fit to analytic models using the XSPEC v. 150 fit- 

1 See jhttp:/ /xmm.esa.int/sas| 

2 http: / /heasarc. gsfc.nasa.gov / docs / xanadu/xspec / index.html 



ting package (|Arnaudl ll996T ). EPIC-pn, MOS1 and MOS2 
spectra were fitted simultaneously, but with normalizations 
varying independently. For timing analysis we used standard 
XANADU/XR ONOS v M tasks and Lomb-Sc argle period 
search method (|Scargle!ll98"2l iPress et al.ll2007l ). 

The data of Chandra observations was processed using 
the CIAO v4.43 threads. We performed standard screening 
of the Chandra data to exclude time intervals with high 
background levels. For each observation, we generated X- 
ray images in the 0.3-7 keV energy band, and used CIAO 
wavelet detection routine wavdetect to detect point sources. 

To estimate upper limits on the quiescent source lumi- 
nosities, the Chandra/ ACIS and XMM-Newton/EPIC count 
rates were converted into energy fluxes in the 0.3-10 keV 
energy range using Web PIMM93, assuming an absorbed 
power law spectral shape with photon index F = 1.2 and 
foreground absorbing column Nn=l.5 x 10 20 cm -2 . 

In the fol lowing analysis we as sume a distance of 817 
kpc for M33 |Freedman et al.ll200ff h All parameter errors 
quoted are 68% (la) confidence limits. 



3 RESULTS 

3.1 Source detection and optical counterparts 

A new X-ray source XMMU J013359. 5+303634 has been dis- 
covered in the data of the 2010 July 9-10 and 11-12 XMM- 
Newton observations of the M33 field (Table [TJ. The es- 
timated source luminosity was ~ 1.4 x 10 37 ergs s -1 , as- 
suming the distance of 817 kpc. XMMU J013359. 5+303634 
was detected during the subsequent, 2010 August 12-13 
XMM-Newton observation, but at much lower luminosity 
level of ~ 5.4 x 10 35 ergs s _1 . We measure the position 
of XMMU J013359.5+303634 to be a = 0l' l 33 m 59.52 s , S = 
+30°36'34.5" (J2000 equinox) with an uncertainty of ~ 1.0" 
(Fig. [IJ. The projected galactocentric distance of XMMU 
J013359. 5+303634 is ~ 3.6' or ~ 850 pc at 817 kpc. The 
analysis of other archival observations of the same field with 
Chandra did not yield source detection with an upper limit 
(2cr) ranging from ~ 5.9 x 10 34 to ~ 5.7 x 10 35 ergs s -1 (or ~ 
24-230 times lower than the observed outburst luminosity), 
depending on the duration of the observation and instru- 
ment used (Fig. [2]). In addition, the source was not detected 
in the 2001-2003 XM M-Newton survey observations of M33 
l|Pietsch et alj|2004h down to an average level of ~ 10 ergs 

3 http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xronos/xronos. html 

4 http://asc.harvard.edu/ciao/ 

5 http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 
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Table 2. Objects from Local Group Survey Galaxies Survey (LGGS) catalog (jM assev et al,lEo06h located within the 3<r error circle of 
XMMU J013359.5+303634. 



LGGS ID 


RA (J2000) 


Dec (J2000) 


V mag 


(B-V) 


(U-B) 


A 




(h:m:s) 


(d:m:s) 








(") 


105063 


01:33:59.47 


+30:36:33.0 


22.40 


0.30 


-0.01 


1.6 


105149 


01:33:59.52 


+30:36:36.0 


20.92 


-0.11 


-0.86 


1.5 


105311 


01:33:59.60 


+30:36:34.3 


21.17 


0.39 


-0.96 


1.1 



The search for the optical counterparts using the images 
from the Digitized Sky Survey (red and blue) did not yield 
any stellar-like object brighter than ~ 20 mag within the 3<r 
error circle of XMMU J013359.5+303634, suggesting that it 
is not a foreground object and probably belongs to M33. 

A cross-correlation with optical catalog based on the 
deep, high resolution Local Group Galaxies Survey (LGGS) 
jMassev et all 120061 ) yielded three possible counterparts of 
~21-22 mag in V band within 3" distance from the X-ray 
position (Table [2 1. Assuming the M 33 optical color excess 
E(B - VP =0.12 H|Massev et all 120071 '). the brightest optical 
object (LGGS 105149) has intrinsic colors (B - V)o « -0.2 
and (U - B) w -0.9 and magnitude Mvo ~ -4 at 817 
kpc consistent w ith being an early B V star jMassev et al.l 
l200fj ; ICoxH2000h . It also falls into Be star domain in the 
(U—B)/(B — V) dia gram, as does its reddening-free Johnson 
Q-ind ex of ~ -0.8 (|Feinstein fc Marracdll979l ; lMassev et all 
l2007t ). The second brightest optical object (LGGS 105311) is 
probably a blend of unresolved two or more stars of different 
spectral class, having (U — B)q ~ —1.0 typical of O to B star 
and (B — V)o ~ 0.3 of late A to early F spectral class jCoxl 
120001) ■ 



ing the periods determined from epoch folding analysis. The 
resulting folded light curves of XMMU J013359. 5+303634 
in the 0.3-10 keV energy band during first two observa- 
tions are shown in Fig. [3] (lower panels). The source demon- 
strates quasi-sinusoidal pulse profiles with additional signif- 
icant peak around phase 0.3 (Fig. [3]). The pulsed fraction, 
defined as (J max - imin)/(-fmax + Imin), where I max and J min 
represent source intensities at the maximum and minimum 
of the pulse profile excluding background photons, was rel- 
atively stable during both Jul. 9-10 and 11-12 observations 
with respective amplitudes of 46±5 and 47±3 per cent (Ta- 
ble(3|. We did not detect any significant change of the pulsa- 
tion period between the 2010 Jul. 9-10 and Jul. 11-12 XMM- 
Newton observations. 

To investigate energy dependence of the source pulse 
profile, we extracted light curves in the soft (0.3-2 keV) and 
hard (2-10 keV) bands for 2010 Jul. 11-12 observation, and 
folded them at the corresponding best pulsation period (Fig. 
3J. Because of the relatively poor statistics, we could not de- 
tect any significant difference between source pulse profiles 
at low and high energies, which have background corrected 
pulsed fractions of 46+4 (0.3-2 keV) and 47+4 per cent (2-10 
keV). 



3.2 X-ray pulsations 

We performed timing analysis of XMMU J013359. 5+303634 
using the data from XMM-Newton/EPIC detectors in the 
0.3-10 keV energy band. After a barycentric correction of 
the photon arrival times in the original event lists, we per- 
formed a Fast Fourier Transform (FFT) and Lomb-Scargle 
periodogram analysis, in order to search for coherent period- 
icities. For the analysis of XMM-Newton data, we used com- 
bined synchronized EPIC-pn and MOS light curves with 2.6 
s time bins to improve sensitivity. We found strong peaks in 
the Fourier spectra and periodograms of data from the 2010 
Jul. 9-10 and Jul. 11-12 XMM observations at the frequency 
of ~ 3.504 x 10 -3 Hz (Fig. upper panels). The strengths 
of the peaks in the individual Fourier spectra (Fig. [3]) cor- 
respond to the period detection confidence of ~ 7 x 10~ 6 
and ~ 3 x 10~ 23 for the Jul. 9-10 and Jul. 11-12 XMM- 
Newton observations jVaughan et al.fl994l ) . The source pho- 
ton count during the 2010 Aug. 12-13 observation was too 
low to perform timing analysis. 

To refine the estimate of the pulsation period, we used 
an epoch folding technique, assuming no period change dur- 
ing individual observations. The most likely values of the 
pulsation period (Table [3]) were obtained fitting the peaks 
in the x 2 versus trial period distribution with a combination 
of Gaussian function and a constant. The period errors in 
Table [3] were computed following the procedure described in 
l|Leahvill987h . Then the source light curves were folded us- 



3.3 X-ray spectra 

The pulse phase averaged XMM- Newton /EPIC spectra of 
XMMU J013359. 5+303634 can be adequately fit with the 
absorbed simple power law model with photon index, T ~ 
1.2 and an equivalent hydrogen column density Nh ~ 1.5 x 
10 21 cm' 2 . The corresponding absorbed luminosity of the 
source in the 0.3-10 keV band is ~ 1.4 x 10 37 ergs s _1 , as- 
suming the distance of 817 kpc. The best-fit spectral model 
parameters of the source are given in Table [3] The measured 
absorbing column A^h is ~2.5 times higher than the Galac- 
tic hydrogen column in th e direction of M33, 6xl0 20 cm -2 
jDickev fc LockmarJlT990T ), consistent with an additional in- 
trinsic absorption within the system and inside the disk of 
M33. 



4 DISCUSSION AND CONCLUSIONS 

The absence of the bright optical counterpart to XMMU 
J013359. 5+303634, its overall X-ray properties and posi- 
tional coincidence with one of the major spiral arms in M 
33 (Fig. [T]) allow us to conclude that it should be located 
outside of our Galaxy and most probably belongs to M 33. 
The X-ray pulsations and spectra of the source are consis- 
tent with that of the a ccreting highly magnetized neutron 
star in a binary system (|White. Swank fc Ho lt 1983; Nagasd 
Il989t iBildsten et alJll997h ■ 
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Table 3. X-ray pulsation parameters and spectral fit information for XMMU J013359. 5+303634. 





Timing Parameters 






POWERLAW*WABS Spectral Model Parameters 




Date, UT 


Period PFo.3-iofeeV' 
(s) (per cent)" 


N H 

(xlO 20 cm~ 2 ) 


Photon 
Index 


Flux 6 


Flux c 


r d 

^X 




x 2 

(d.o.f.) 


Jul. 9-10 
Jul. 11-12 


285.427 ± 0.089 
285.365 ± 0.047 


46 ±5 

47 ±3 


14 ±5 

15 ±2 


i 99 +0.13 
1 - zz -0.12 

1.17 ±0.06 


1.60 ±0.07 
1.69 ±0.04 


1.76 ±0.14 

1 ocr+0.09 
1 -°°-0.08 


1.28 
1.35 


1.41 
1.48 


67.6(69) 
190.2(143) 



a - pulsed fraction in the 0.3 — 10 keV energy band 
b - absorbed model flux in the 0.3 — 10 keV energy range in units of 10 — 13 erg s _1 cm -2 
c - unabsorbed model flux in the 0.3 — 10 keV energy range in units of 10 — 13 erg s — 1 cm -2 
d - absorbed luminosity in the 0.3 — 10 keV energy range in units of 10 37 erg s _1 , assuming the distance of 817 kpc 
e — unabsorbed luminosity in the 0.3 — 10 keV energy range in units of 10 37 erg s _1 , assuming the distance of 817 kpc 



In the high-mass system scenario, the 285.4 s pulse pe- 
riod places XMMU J013359.5+303634 squarely among Be 
systems or underfilled Roche Lobe supergiant systems o n 
the Corbet diagram (|Corbetl [l98r3 ; ICharles fc Cod l200fj ). 
The transient behavior of the source favors the Be X-ray 
binary interpretation, in which the X-ray outbursts are at- 
tributed to either periastron passages of the compact com- 
panion (Type I outbursts) or changes in the Be star circum- 
stellar disk that might occur at any binary phase (T ype II 
outbursts) jStella et alJll986t iNegueruela et al.ll200ll ). Pos- 
sible optical counterparts to XMMU J013359. 5+303634 are 
all fainter than 20 mag in V band (Table [2]), with absolute 
V magnitudes between ~ —2 and ~ —4 at 817 kpc, lend- 
ing further support to a non-supergiant system interpreta- 
tion. In addition, the brightest optical object located within 
the error circle of XMMU J013359. 5+303634 has properties 
consistent with being an early B V star and possibly be- 
longing to Be class. Assuming the distance of 817 kpc, the 
0.3-10 keV X-ray luminosity of XMMU J013359. 5+303634 
measured in Jul. 9-12 observations (-L0.3-10 ~ 1.3 x 10 37 
ergs s~ x ) also f alls within the range observed for Be X-ray 
binary pulsars (|Bildsten et al.lll997l ). 

If XMMU J013359.5+303634 is indeed a Be system, its 
binary period, P or b estimated from the Corbet diagram falls 
in the ~ 100 — 200 day range. The observed time span be- 
tween XMM-Newton observations when the source was de- 
tected allows us to put an additional limit on the orbital 
period. Assuming that we observed a single Type I outburst 
lasting less than one orbital cycle requires P or b > 35 days. 

Another possibility is that XMM U J013359. 5+30363 4 
is a transient low- mass binary system (|Bildsten et al.|[l997l ). 
In that scenario, the transient outbursts could be explained 
as a result of viscous-thermal instability of the quies- 
cent accretion disk. The relatively high pulsation period 
of the source probably excludes short orbital period sys- 
tems, but the possibility of the very long period system with 
red giant companion similar to Galactic source GX 1+4 
ijChakrabartv fc Rochel Il997h can not be completely ruled 
out. 

The proposed high-mass nature of the system implies 
its relatively young age, consistent with its location within 
one of the spiral arms of M 33 (Fig. [T]). Since the Be in- 
terpretation still remains preliminary, optical identification 
is essential to determine the nature of the system. To fa- 
cilitate it, follow-up monitoring observations with Chandra 
and XMM-Newton are needed to improve source localiza- 
tion, study long-term evolution of its X-ray properties and 



possibly putting constraints on orbital parameters of the bi- 
nary system. 

The number of X-ray pulsar candidates located out- 
side our own Galaxy and Magellanic Clouds is still very 
small with only four systems known so far, including XMMU 
J013359. 5+303634. All of them were discovered as tran- 
sient/recurrent sources with X-ray luminosities in the ~ 
10 37 — 10 39 ergs s _1 range, and probably belong to a lu- 
minous Be/X-ray binaries class. The discovery of XMMU 
J013359. 5+303634 is yet another demonstration that deep 
monitoring observations with XMM-Newton and Chandra 
can detect these and similar pulsating sources (both tran- 
sient and persistent) in the galaxies located up to a few Mpc 
away. Since the nearby normal galaxies with recent and on- 
going star formation provide the environment where X-ray 
pulsars are expected to be plentiful, careful analysis of the 
existing archival and new observations of these galaxies have 
a potential to significantly increase the statistics of X-ray 
pulsar systems and improve our understanding of their na- 
ture and connection to the underlying stellar populations. 
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Figure 1. Upper panel: Combined 0.3-7 keV XMM/EPIC-MOS 
image covering a 12' X 12' central region of M33, taken on Jun. 
11-12, 2010. The position of new pulsar XMMU J013359.5-303634 
is marked with an arrow, and the position of M33 nucleus is 
shown with a white cross. Lower panel: Optical (B-band) im- 
age of t he central region of M33 from the Local Group Galaxies 
survey dMassev et alj|200fj|) . The region covered by X-ray image 
in the upper panel is shown with black rectangle. The position of 
XMMU J013359. 5-303634 is shown with black circle of 5" radius. 
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Figure 2. Upper panel: X-ray light curve of XMMU 
J013359. 5+303634 in the 0.3-10 keV energy band based on the re- 
sults of 2010 XMM-Newton observations. Lower panel: Long-term 
light curve of XMMU J013359.5+303634 based on the results of 
2000-2010 Chandra and XMM-Newton observations. The source 
luminosities were calculated assuming the distance of 817 kpc. 
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Figure 3. Upper panels: Power spectra of XMMU 
J013359. 5+303634 obtained using the data of 2010 Jul. 9- 
10 and 11-12 XMM-Newton/EPIC (EPIC-pn, MOS1 and MOS2 
detectors combined) observations in the 0.3-10 keV energy band. 
Lower panels: Corresponding background corrected source X-ray 
pulse profiles folded with most likely pulsation periods (Table ■ 



XMM-Newton discovery of transient X-ray pulsar in M33 7 



XMMU J013359.5 + 303634; Jul. 11 — 1 3. 2010; XMM/EPIC 
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Figure 4. Background-corrected normalized X-ray light curves of 
XMMU J013359.5+303634 during the 2010 Jul. 11-12 observation 
folded at the corresponding best period (Table[3]l in the 0.3-2 and 
2-10 keV energy bands (upper and middle panels) along with their 
hardness ratio (bottom panel). 
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Figure 5. EPIC count spectra and model ratios of XMMU 
J013359.5+303634 during the 2010 Jul. 11-12 observation. The 
EPIC-pn data is plotted with filled circles, while EPIC-MOS1 
and MOS2 data is shown with filled and open squares respec- 
tively. The best-fit absorbed power law model approximation of 
EPIC-pn, MOS1 and MOS2 data is shown with solid, dotted and 
dashed histograms. 



